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Abstract. For partially wetting, ellipsoidal colloids trapped at a fluid interface, their effective, interface- 
mediated interactions of capillary and fluctuation-induced type are analyzed. For contact angles different 
from 90°, static interface deformations arise which lead to anisotropic capillary forces that are substantial 
already for micrometer-sized particles. The capillary problem is solved using an efficient perturbative 
treatment which allows a fast determination of the capillary interaction for all distances between and 
orientations of two particles. Besides static capillary forces, fluctuation-induced forces caused by thermally 
excited capillary waves arise at fluid interfaces. For the specific choice of a spatially fixed three-phase 
contact line, the asymptotic behavior of the fluctuation-induced force is determined analytically for both 
the close-distance and the long-distance regime and compared to numerical solutions. 
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phenomena, random processes, noise, and Brownian motion 



05.40.- 



Fluctuation 



in 

00 

O 
00 
O 



X 



1 Introduction 

Colloidal particles trapped at fluid interfaces exhibit effec- 
tive interactions which broadly can be classified into "di- 
rect" interactions (e.g., of electrostatic or van-der-Waals 
type) which also exist for colloids in bulk solvents but 
they are modified at interfaces. Additionally, the presence 
of an interface gives rise to interactions mediated by de- 
formations of that interface, thus these are absent in bulk 
solutions of colloids. One speaks of capillary interactions 
[1], if these deformations are static, and of fluctuation- 
induced interactions if the deformations are caused by 
thermal fluctuations. For spherical colloids at free inter- 
faces, the formation of self-organized structures is mainly 
governed by the direct interactions whereas for particles 
of nonspherical shape capillary interactions appear to be 
dominant [2]. 

Static interface deformations arise if the colloids expe- 
rience forces in the direction parallel to the interface nor- 
mal (e.g. if they are pushed into the lower phase) and/or 
stress distributions act on the interface [3]. For microcol- 
loids of sizes less than 10 yum, the omnipresent gravita- 
tional force on the colloids can be neglected. However, 
meniscus deformations also arise in conjunction with di- 
rect interactions (like electrostatic forces) which lead to 
forces and stresses on colloids and interface, respectively. If 
the system "colloids -t- interface" is mechanically isolated 
(usually this applies for colloid experiments in a Langmuir 
trough or on large droplets), then the force on the colloids 
directed vertical to the interface is balanced by the total 



force on the interface (obtained by integrating the stress 
distribution over the interface area) [HEUS]- The ensuing 
capillary interactions decay with a power-law in the inter- 
colloidal distance d, e.g., in the case of charged colloidal 
spheres at air-water or oil-water interfaces they are at- 
tractive, oc fi~^, but for large d they are usually weaker 
than the direct electrostatic repulsion which is also oc 
[7ll8l[9]. - On the other hand, in the absence of forces on 
the colloids and stresses on the interface, static interface 
deformations can also be induced by an anisotropic colloid 
shape; more precisely if the colloid is not symmetric with 
respect to rotations around any axis through the colloid 
which is parallel to the normal on the undisturbed inter- 
face. Young's equation requires that at the three-phase 
contact line the angle between the local interface nor- 
mal and the local normal on the colloid surface is given 
by the contact angle 6. Thus for an anisotropic colloid 
this condition cannot be met if the interface remains flat; 
the contact line will not be located in the plane of the 
undisturbed interface. The associated interface deforma- 
tions around one such colloid can be calculated in terms of 
a two-dimensional multipolar expansion [lOj . For asymp- 
totically large distances from the colloid, the leading non- 
vanishing multipole is in general the quadrupole, since 
monopole and dipole are absent through the conditions 
of force and torque balance. The interaction energy be- 
tween two quadrupoles depends on the colloid orientation 
in the interface plane and decays according to a power law, 
oc c?~^ . Experimentally, it has become possible to produce 
anisotropic microcolloids of controlled shape for investiga- 
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tions of their self-assembly at fluid interfaces: anorganic 
nanorods [11] , bended disks [12] or rotationally symmetric 
ellipsoids (spheroids) with aspect ratios up to 10 |T3l[T4] . 
For the ellipsoids, results on the effective pair potential 
for intermediate distances d |13| indicate a strong orienta- 
tion dependence not captured by the aformentioned lead- 
ing quadrupolc interaction. Ellipsometric measurements 
of the interface deformation around one particle are con- 
sistent with a quadrupolar pattern [15] which, however, 
appears to be deformed considerably for stretched ellip- 
soids. 

The calculation of the full interface profile around trap- 
ped ellipsoids and the evaluation of the associated capil- 
lary interaction energy can only be done numerically, save 
for small eccentricities e where an expansion in terms of 
e is possible |16| . Similarly, an evaluation of the inter- 
face profile is simplified if for colloidal spheres only small, 
roughness-induced deviations from the circular shape of 
the contact line are assumed [TTlfTS]. For ellipsoids with 
eccentricities not close to zero, however, the position of the 
contact line is not given a priori but has to be determined 
self-consistcntly through Young's equation. This task is a 
variant of a free boundary problem which in general poses 
problems in terms of speed and efficiency of a numerical 
algorithm aimed at solving it. Below, we show that a per- 
turbative solution can be attained through expansion of 
an appropriate free energy functional and subsequent min- 
imization which leads to a problem with fixed boundaries 
(Subsec. 12. ip . The calculation of the capillary interaction 
between two ellipsoids is thus greatly simplified and re- 
sults can be obtained for the full range of distances and 
orientations in a configuration with two ellipsoids (Sub- 
sec. E2]). 

The capillary interactions between two ellipsoids are 
absent if the contact angle 6 equals tt/2, and they be- 
come significantly smaller for ellipsoid sizes approaching 
the molecular length scale. In these circumstances, fluc- 
tuations around the equilibrium interface position are ex- 
pected to influence the effective interactions noticeably. 
Since the interface fluctuations are of thermal nature, the 
scale of the ensuing interactions will be given by fceT, 
the thermal energy. Within a coarse-grained picture, the 
properties of fluid interfaces are very well described by an 
effective capillary wave Hamiltonian which governs both 
the equilibrium interface configuration and the thermal 
fluctuations (capillary waves) around this equilibrium (or 
mean-field) position. As postulated by the Goldstone the- 
orem the capillary waves are long-range correlated. The 
interface breaks the continuous translational symmetry of 
the system, and in the limit of vanishing external fields 
- like gravity - it has to be accompanied by easily ex- 
citable long wavelength (Goldstone) modes - precisely the 
capillary waves. The fluctuation spectrum of the capillary 
waves will be modified by colloids trapped at the interface 
and therefore leads to fluctuation-induced forces between 
them. In that respect, colloids at fluids interfaces appear 
to be a possible realization of a two-dimensional system 
exhibiting the Casimir effect. For spheres and disks, these 
forces have been calculated in Refs. [T9ll20j. The large- 



distance behavior of these forces depends sensitively on 
the boundary conditions at the three-phase contact line 
whereas for close distances a strong attraction similar to 
van-der-Waals forces has been found, independent of the 
type of boundary condition. In the case of colloidal rods 
the asymptotic behavior of the fluctuation-induced force 
has been evaluated in Ref. [21] and shown to lead to an 
oricntational dependence. Furthermore we note that there 
is numerous work on the force between inclusions on mem- 
branes where the membrane shape fluctuations take the 
role of capillary waves, see, e.g., Refs. [^b22j . 

In Sec. [3] below, we consider the specific case of two el- 
lipsoids trapped at an interface with a pinned contact line. 
The fluctuation-induced force between two such ellipsoids 
corresponds to the Casimir force with Dirichlct boundary 
conditions. We present numerical results for the whole dis- 
tance regime for selected orientations and compare to the 
leading terms in analytic expansions valid for the close- 
distance and the long-distance regime. 

Finally, Sec. [4] contains a discussion of the results. 

2 Static interface deformations and capillary 
interactions 

We consider cigar-shaped ellipsoids of micrometer size 
with half-axes (a, 6, 6) and a > b {e — {1 ~ &^/a^)^/^ is 
the eccentricity) trapped at an air-water interface with 
surface tension 7. The surface tensions of the ellipsoid 
with air and water are denoted by 71 and 711, respec- 
tively. The contact angle (Young's angle) is defined by 
cos 9 = (71 — 711) /j. The free energy of the system shall be 
given only by the surface free energies of the three involved 
interfaces (ellipsoid-air, ellipsoid-water and air-water), 
i.e., we neglect gravitational effects which are negligible 
for micrometer-sized particles and also possible electro- 
static effects which arise through the ubiquituous surface 
charges on real colloids. (We will comment upon electro- 
static effects in Subsec. [2^ . 

2.1 Capillary deformation around a single ellipsoid 

For 6 in the range between and 180° and in the absence 
of line tensions, the most stable configuration is given by 
the ellipsoid positioned flat on the interface (Fig. [T]) , be- 
cause in this configuration the amount of displaced area of 
the air-water interface is maximal^ The meniscus defor- 
mation with respect to the plane z = is denoted by u{r) 
where r is a two-dimensional vector in the plane z = 
and the vertical position of the ellipsoid center is given by 
h. In order to determine the equilibrium deformation and 
vertical ellipsoid position, we seek an expansion of the free 
energy around a reference configuration characterized by 
the deformation Uiof (r) and the vertical position hj-ot such 
that 

U = Urcf + w , h = hici + . (1) 

^ This can be shown by using the simplifying assumption 
that the interface around the ellipsoid remains flat |23j . 
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Fig. 1. (a) Deformation of a fluid interface around an ellipsoid 
(satisfying Young's equation locally), with Aumux denoting the 
maximal difference in meniscus elevation along the three-phase 
contact line. The part of the ellipsoid protruding out of the wa- 
ter is shown unshaded and the part of the ellipsoid immersed in 
water is depicted in darker shade, (b) Side view of the ellipsoid 
with the geometrical definition of the contact angle 9. 

The properties of the reference configuration will be de- 
termined below. The change in free energy with respect 
to this reference configuration is given by 



- jiAAi + -/uAAu 



(2) 



where Z\Amen denotes the change in meniscus area and 
ZiAi[ii] denotes the change in contact area between the 
ellipsoid and air or water, respectively. We split this free 
energy difference into two parts: 



AJ^ = T,r 



(3) 
(4) 



where the meniscus free energy .Fmen denotes the differ- 
ence of the air-water interfacial energy integrated over 
that part of the plane z = which is given by the area of 
the meniscus Wref of the reference configuration projected 
onto it (denoted by S^ci)- The boundary free energy .Fb 
includes all remaining terms. We denote by S the area 
of the meniscus u projected onto the plane z = and 
Zou(a;,J/;/i) describes the surface equation of the ellipsoid 
which depends on its vertical position /i as a parameter. 
With these definitions 



7/ (i^r cos 6*^1 + (Vzcii) 



(5) 



At this point we perform a Taylor expansion up to second 
order in the meniscus deformation for the free energy con- 
tribution JFinen and separately. For .Fmen, this is equiv- 
alent to the small-gradient expansion |Vw|,|Vurof| ■C 1 
and results in 



.f + v]f-(yu,,tf) . (6) 



The expansion of the boundary part J^i^ is somewhat more 
involved. Since the area of the domain S'ref\S' will be of 
second order in the meniscus deformation, we can approx- 
imate \/l + (V[Mrcf + v])^ w 1 in Eq. ^ and find 



2tt 



dr r 

ro.iaf 



(cos^yTTTv^ 



(7) 



Here the functions ro,rof(0) and r((^) = ?'o.ref + Ar{(f)) 
parametrize the polar radius of the boundaries dS^et and 
dS, respectively, i.e., they correspond to the projected 
three-phase contact lines formed by the reference menis- 
cus Urcf and the arbitrary meniscus u. Subsequently we 
perform a functional expansion of Eq. ([7]) with respect to 
displacements of the contact line position = v{r {(/))) — 
Ah which are constrained to lie on the ellipsoid surface: 



Svd, 



(8) 



■e=o 



2 J J Sv^Sv^' 



v=0 



The position of the reference contact line is fixed by the re- 
quirement that the reference configuration minimizes Ji,, 
i.e., 



dr{(j)) 5!F\y 







r(0)=ro,rof (0) 



5=0 ^^i't)) 

which leads to the following condition on ro.rcf(0): 
cos^? — 



y/l + [VZell(ro,rcf (</>))] ^ 



(9) 



(10) 



Geometrically, Eq. (|T0|) expresses the condition that at the 
reference contact line the angle between the unit vector in 
z-direction and the ellipsoid normal is given by 9 which is 
a reasonable first "guess" of the equilibrium contact line 
position. The solution to this equation yields ellipses for 
the projection of the reference contact line with half axes 
h' = b\ sine\ and a' = b'/{l - e^[l + cos^ ^(l - e^)])!/^. 
In order to fully specify the reference configuration, we 
require that the reference meniscus is a surface with mini- 
mal area, i.e, to first order it fulfills Auj-of = (consistent 
with Eq. ([6])) with the pinning condition Wref |ro „t(<7i) = 

^ell(ro,ref((^); /iref)- 

The second term on the r.h.s. of Eq. ([8|) gives us the 
approximation for the boundary free energy used in the 
following: 



R{cp) 



7 



d<pR{cj>)il 
sin^ 9 



?'0,rcf (0)^(1 - C0S2 



(11) 

(12) 



Since it is a local functional on the projected reference 
contact line, this boundary free energy can be viewed as 
the total free energy cost in shifting the contact line with 
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respect to the reference state. Thus we see that the min- 
imization of the free energy AJ- = J-mcn + J'h, Taylor 
expanded to second order with JFmon given by Eq. ([6]) 
and Th by Eq. (fTTj) . with respect to v leads to the lin- 
earized Young-Laplace equation Aw = (which also im- 
plies Au = through Aiticf = 0) with the local boundary 
condition on the projected reference contact line ro,rcf('?^') 



d{Ure{ + v) 

dn 



Ah) . (13) 



Here, d/dn denotes the outward normal derivative on dS^^i 
^o,ref('/') and M is the line differential on dS^ef- Through 
minimization of AJ- with respect to Ah^ the vertical po- 
sition of the colloid is fixed by the condition 



Ah 



d4)R{(t)) w(ro,rof (0)) 



d(l)R{(j)) 



(14) 
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Fig. 2. Comparison of the elliptic and polar multipole ex- 
pansion of the equilibrium meniscus around an ellipsoid with 
aspect ratio a/b — 5 and contact angle 6 — 66°. The elliptic 
multipole series is defined in Eq. (|17p and the polar multipole 
series (where we have used ro — a) is given by Eq. (|18|) . Note 
that because of symmetriy reasons all sine moments Bm and 
BE, are zero. 



This condition implies that the solution does not depend 
on the choice of /irof , the vertical position of the ellipsoid 
in the reference configuration. An arbitrary shift in h^^f 
will be compensated by a corresponding negative shift in 
Ah, as can be shown from Eqs. and p4|) . 

At this point we want to remark that the reduction 
of the original capillary problem (with unknown contact 
line) to the solution of the Laplace equation with a lo- 
cal boundary condition at a fixed boundary (the projec- 
tion of the reference contact line) is a great simplification 
and speeds up numerical solutions enormously. Further- 
more we note that the technique of splitting the free en- 
ergy into a meniscus and a boundary part with subsequent 
Taylor expansion has already been introduced in Rcf. [24| 
for the problem of capillary deformations around colloidal 
spheres. 

Since the boundary curve ro,rcf(0) is itself an ellipse, 
a numerical solution for the meniscus deformation u is 
most conveniently performed using elliptic coordinates s, t 
whose relation to cartesian coordinates is given by x = 
a cosh s cost and y ~ a sinhssint. Isolines of constant s 
are ellipses, the condition that the boundary curve ro^rcf (</*) 
is such an isolinc at s = sq leads to the relations a = a'e' 
and So = acosh(l/e') where e' is the eccentricity of the 
elliptic boundary curve. The Laplace equation in elliptic 
coordinates is given by 

Au(s, t) = — i^— + —j u{s, i) = , (15) 

H = aV sinh^ s + sin^ t , (16) 

and its solution is given by the expansion 

u{s,t) = Ao — + (17) 
So 

e-ni{s-so) cos{mt) + Bm sm{mt)] , 

m>0 



where A,„ and Bm denote elliptic multipole moments of 
order m. Comparing to the general solution in polar coor- 
dinates, 

u{r, (j)) = AP\n—+ (18) 
^ [A^m cos(™0) +i?P sin(™0)] , 

m>0 

we note that an elliptic multipole of order m is a super- 
position of polar multipoles of order n > m. In terms of 
the expansion given in Eq. ([T7)) . the problem reduces to a 
set of coupled linear equations for the multipole moments. 
It is usually sufficient to take into account multipoles of 
order m < 50 for a very precise solution. 

In Fig.[2]we compare the convergence of the elliptic and 
the polar multipole series for the meniscus deformation 
around an ellipsoid with aspect ratio a/6 = 5 and contact 
angle 6 — 66°. The leading multipole is the quadrupole 
(?n 2). The vanishing monopole is related to the fact 
that no contact line force acts on the ellipsoid which we 
have ensured by minimizing the free energy with respect to 
the vertical position of the ellipsoid. The dipole moments 
are also zero since in equilibrium there is no torque acting 
on the ellipsoid. While the elliptic multipole expansion 
converges sufficiently fast for all practical purposes (e.g., 
A30/A2 ~ 10~^), the polar multipole series is very badly 
convergent (e.g., ^gp/Aj ~ 0.03). This is due to the fairly 
large aspect ratio; for a/b ^ 1 the polar multipole series 
yields rapid convergence |16| . 

Due to the approximate forms of the meniscus and 
boundary free energy functionals (Eqs. © and (jlip). Young's 
condition is fulfilled only approximately for the correspond- 
ing minimum configuration. In general, the approximation 
becomes exact in the limit 9 —>■ 90° (for arbitrary aspect 
ratio a/b) or a/b ^ I (for arbitrary contact angle 6). The 
projected contact line ro{4>) ^ ?'o,icf (</>) of the approximate 
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Fig. 3. The ratio cos Q' j cos Q for the two aspect ratios 2 and 5. 
The angle d' is the contact angle pertaining to the approximate 
solution (varying along the contact line contour, parametrized 
by the elliptic angle i) and Q = 66° is the input contact angle. 
The inset shows the projection of the ellipsoid contour and the 
reference contact line on the interface plane. For some values 
of t, the location of the corresponding point on the contact line 
is given. 

solution can be determined numerically by the intersection 
of the solution (fT7|) with the ellipsoid. The likewise numer- 
ically determined contact angle 6' varies along Tq (</)), for 
an example see Fig. [3] where the deviations from Young's 
law are shown for an input contact angle 9 = 66° and the 
two aspect ratios a/b = 2 and a/b = 5. For the smaller 
value of the aspect ratio, the deviations are small. For 
a/b — 5 however, larger deviations occur which are local- 
ized very close to the tips (in a domain t < 10°, see the 
inset of Fig. [3|) . Closer inspection reveals that at the tips 
and on the contact line ro(0), the small-gradient approx- 
imation (Vit)^ ^ 1 breaks down (however, it still holds 
approximately on the reference contact line ro,rof(0))- 

The meniscus shape strongly depends on the aspect 
ratio and the contact angle. We have investigated the in- 
fluence of these parameters on the maximum height differ- 
ence Ziumax along the contact line, with the results shown 
in Fig. m This provides us with a quick estimate on the 
strength of the interaction between two ellipsoids since in 
linearized theory one can expect that the amplitude of the 
capillary interaction energy is approximately proportional 
to (Z\uniax)^- Depending on the precise value of the aspect 
ratio, Ziitmax attains its maximum for contact angles be- 
tween 40° and 55°, i.e. in the experiments of Ref. [15] . 
where for the used polystyrene ellipsoids a contact angle 
of around 40° was determined, the capillary deformation 
around the particles is close to its maximum as compared 
to particles of the same shape but with different 6. 



2.2 Capillary interaction between two ellipsoids 

For large distances d, the capillary interaction between 
two ellipsoids is determined by the quadrupole. For small 
eccentricities (a ~ & = tq) the interaction energy has been 



Fig. 4. Maximum meniscus height difference Z\Mniax(S) as 
function of the contact angle 9 for different aspect ratios a/b. 
Climax attains its maximum for contact angles between 40° 
and 55°, while it is zero for neutrally {6 = 90°) or completely 
{9 0) wetting particles. 

determined in Rcfs. [TTlfTU] and reads 

U^:f = -3^7(Ziu„,ax)' cos(2cc>i + 2cc>2) , (19) 

where Ui and are the polar orientation angles of the 
ellipsoids in the interface plane with respect to the dis- 
tance vector between their centers (see Fig. [5] (a)). The 
simple energy estimate in Eq. (|19p predicts that the cor- 
responding forces between two ellipsoids approaching each 
other side-by-side (wi = W2 = 90°) or tip-to-tip (wi = 
LU2 ~ 0°) are attractive and equal. However, an experi- 
mental estimate of these interaction forces [l^ revealed 
that the attractive force in the side-by-side configuration 
varied as Fcap oc — d^^-^ whereas for the tip-to-tip con- 
figuration it varied as i^cap cx — in accordance with 
Eq. ()19|1 . The measurements were performed over a lim- 
ited distance range d/a < 4 for aspect ratios a/b ranging 
from 3 to 4.3. 

Thus, for such aspect ratios the quadrupole approxi- 
mation appears to be insufficient. We can apply the for- 
malism developed in the previous subsection also to the 
evaluation of the interface deformation and associated cap- 
illary energy for two ellipsoids. The meniscus part of the 
free energy (Eq. (jj])) remains unchanged save for the ex- 
tend of the integration domain S'lof : here it consists of the 
whole plane with the two ellipses enclosed by the projected 
reference contact line cut out, S'rct = \ Uf^j^S'i^rof The 
boundary free energy for each ellipsoid (Eq. pTjl ) has to 
be amended by additional terms related to the appearance 
of an additional degree of freedom given by the angles (Xi 
of the long ellipsoid axis with the plane z = 0. For fixed 
distance d and orientation angles iO\ and W2, the presence 
of the second ellipsoid causes the first ellipsoid to dip with 
its tip into the water and vice versa, i.e., the free energy 
has to be minimized with respect to rotations around an 

We follow the convention of Ref. [24] and denote with a hat 
all quantities pertaining to the two-ellipsoid configuration. 
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Fig. 5. Top view on two ellipsoids at an interface, (a) The 
polar angles uji and u)2 specify the direction of the long half axis 
of the contact line ellipse with respect to the distance vector 
between the centers of the two ellipses. The polar angle (/f>i on 
the reference ellipse 9Si,rof and the parametrisation (si,ti) — 
(si(s2, ^2), ii(s2, ^2)) of flSi^rof in terms of elliptic coordinates 
with respect to the reference ellipse 9S2,ref is indicated, (b) 
Ri and R2 are the radii of the curvature on those ellipse points 
whose distance is the minimal distance h. 

axis in the interface plane which is perpendicular to the 
symmetry axis. The additional terms in the boundary free 
energy can be determined as before in Sec. l2.1l bv perform- 
ing a Taylor expansion in terms of both, displacements 
of the contact Une height Vi^^- = v{ro{4>i)) — Ahi and 
changes of the tilt angle ai around the reference configu- 
ration Wi^^. = and ai — (i.e., in the reference configu- 
ration the ellipsoids are positioned flat on the undisturbed 
interface). This procedure results in the expression 




+ 2e^xR,a{(l){) [v - Ah,] a, + Rc,c.{<j),) a^,} (20) 

for the total boundary free energy of two ellipsoids. The 
lengthy expressions for the coefficients Rzz , Rza and Raa 
are given in the appendix, see Eqs. ([37]) - ((39)) . (Note that 
i?,,(0) = Ri^), Eq. m-) 

For two ellipsoids, the meniscus free energy of the ref- 
erence configuration, 

-^mcn,rcf(rf,Wl,t^2) ~ / d^a; (VUrof )^ , (21) 

^ "'Srot 

obviously depends on the configuration variables {d, wi , 0^2} 
The reference meniscus fulfills Auref = and the bound- 
ary condition Urof|as,.„f = 2Eii(ro(0i); /li.ref )■ To include 
all dependence on the configuration variables into the to- 
tal free energy, it is suitable to consider the free energy 



difference with respect to the reference state with d 00: 

AT = f^cxv{d,uji,U2) - ^men,iei{d ^ cx),a;i,a;2) + 

Tb{d,uji,uj2) (22) 

d'^x{Vu)'^ +Th{d,uJi,uj2)+con.st. (23) 

Consequently, the capillary potential L^cap can be defined 
as the free energy difference 

U,,p = AT{d,uji,Lj2) - AT{d ^ ^,iJi,Lj2) . (24) 

The equilibrium configuration minimizes the free en- 
ergy in Eq. ([23]) and can be calculated similarly as in 
Sec. 12.11 Minimizing Th with respect to Ahi and ai pro- 
vides the equilibrium height and orientation of ellipsoid 
i and leads to expressions analogous to Eq. (|14p (since 
terms coupling Ahi and a,; vanish). In contrast to the 
single colloid case, the equilibrium meniscus w(x) for two 
ellipsoids has to fulfill boundary conditions at both refer- 
ence ellipses ■ They arise by minimizing the bound- 
ary free energy in Eq. ([20| and contain an additional, 
ai-dependent term as compared to the boundary condi- 
tion for the single ellipsoid (Eq. p3)) ). For the numerical 
determination of the equilibrium meniscus profile, the su- 
perposition ansatz u(x) = ui{si,ti) -1-^2(52,^2) is used. 
Thereby, the functions Ui are given by the expansions of 
the meniscus around a single colloid into elliptic multi- 
poles (Eq. pT)) ). and {si,ti) are elliptic coordinates with 
respect to the center of the reference ellipse Si^ef (see 
Fig. [5|). Through the boundary conditions a set of cou- 
pled linear equations for the multipole moments is derived 
which can be solved with standard numerical methods. 

As an example for the results, the capillary force i^cap = 
-~dUcap/dd in direction of the distance vector between the 
ellipsoids is shown in Fig. [51 for an aspect ratio a/b = 5 
and a contact angle 9 = 66°. The force considerably de- 
viates from the quadrupole form for d/a < 4, i.e., in the 
region where the experimental measurements of Ref. |13| 
have been performed. For these distances, the force does 
not follow a power law but a fit to an effective power-law 
would clearly yield an exponent > — 5 in the side-by-side 
configuration, but also an effective exponent < —5 for 
the tip-to-tip configuration. The reason for this behav- 
ior appears to be that the capillary deformation around 
one ellipsoid is dominated by the elliptic quadrupole, i.e., 
closer to the ellipsoid the capillary deformation has sub- 
stantial contributions from polar multipoles higher than 
the quadrupole (see Fig. [2|) which also influence the pair 
interaction considerably. 

For the parameters used in Fig. [5] the asymptotic cap- 
illary potential, given by the quadrupole form U^^p'^ = 
-Uo (a/d)'^ cos(2wi -I- 2tJ2), one finds the amplitude Uq « 
7 X 10® fceT for ellipsoids with long half axis a = 10 fim at 
the air-water interface. Usually, the experimentally used 
ellipsoids are charge-stabilized which leads to an asymp- 
totically isotropic dipolar repulsion Ud = J7o,ci (a/d)^. Us- 
ing the results of Ref. [9] , one can estimate the amplitude 
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Fig. 6. Capillary force -Fcap = —dUca.p/dd in units of 7a 
between two ellipsoids with aspect ratio a/b = 5 and contact 
angle 9 = 66° approaching each other side-by-side or tip-to- 
tip. 



of the electrostatic repulsions as [/o,ci ~ 10'^ fceT (with a 
charge density of 1 electron per nm^ and ultrapure water). 
Thus the electrostatic repulsions arc completely unim- 
portant compared with the capillary potential; only for 
distances d > 10^ a the directional capillary attractions 
would be overpowered by the electrostatic repulsions. 



3 Capillary waves: fluctuation-induced 
interactions 

As discussed in the introductory section, the scale of the 
fluctuation-induced interaction energies is ~ ksT such 
that they are observable only if static capillary interac- 
tions are (almost) absent. Since for 9 = 90° the capillary 
interactions are identically zero (see the previous section) 
we discuss in this section the exemplary case of two el- 
lipsoids with their centers at z = and their three-phase 
contact lines (ellipses with half axes a and b) also located 
in the plane of the flat interface, corresponding to the equi- 
librium position of an ellipsoid with contact angle 6 = 90°. 
Furthermore we assume that the contact line is pinned 
and the position of the ellipsoid is fixed by some external 
means. This corresponds to a Dirichlet boundary condi- 
tion for the meniscus at the boundary of the projected 
meniscus area S: u{dS) = 0. The free energy cost for 
small-gradient fluctuations of the interface position u(r) 
around its mean u = is given by the capillary wave 
Hamiltonian: 



He 



= 1 1 d^' 



(25) 



which (with Ac — ^ 00) has already been used for evalu- 
ating the free energy of static interface deformations (see 
Eq. ([n])). In Eq. (|25p . Ac denotes the capillary length which 
is usually much larger than the extensions of microcol- 
loids, nevertheless it is necessary to keep the associated 
free energy contribution (stemming from gravity) through- 
out the calculations since it ensures the stability of the 



interface |20j . The capillary wave Hamiltonian depends 
on the position and orientation of the ellipsoids through 
the integration domain S which encompasses the whole 
plane z — Q with the ellipses enclosed by the contact lines 
cut out. Therefore also the the free energy T{d,uji) = 
—k^T lnZ(d, Wi) depends on the distance d and the ori- 
entation angles oji and L02 (see Fig. [5] (a)). The partition 
function Zid^LOi) is obtained by a functional integral over 
all possible interface configurations u] the boundary con- 
dition u\QSi = on the two contact lines dSi and 882 
is included by (5-function constraints, as introduced in 
Ref. [261 : 



Vu cxp < — 



He 



n n '^["(^')]- 

i=i x.eaSi 

(26) 

Zq is a normalization factor such that Z{d — > 00) ~ 1. The 
i5-functions in Eq. (|26[) can be removed by using their in- 



kBT 



tegral representation via auxiliary fields defined on 

the contact lines dSi pslES] . This enables us to integrate 
out the field u leading to 



Z - Zg-i / Y[Vij, exp 



knT 



E 



de,. 



9S, 



dlj X 



i/'»(x,)G(|x, -Xjl)V'j(xj) 



(27) 



Here, we introduced the Green's function 
G(x) = ii:o(|x|/Ac)/(27r) of the operator (- A+A^^) where 
Kq is the modified Bessel function of the second kind. In 
this form, the fluctuation part resembles 2d screened elec- 
trostatics: it is the partition function of a system of fluc- 
tuating charge densities ipi residing on the contact lines. 

In the intermediate asymptotic regime a, 6 <C d <C Ac 
the free energy associated with this partition function can 
be calculated through an expansion of the auxiliary flelds 
ipi in terms of elliptic multipoles and an expansion of 
Green's function in terms of 1/d using elliptic coordinates. 
Details of the lengthy calculations will be reported else- 
where. The resulting free energy can be written as an ex- 
pansion T/ik^T) = fo + !2ld^ + jijd^ + . . . with the two 
leading coefficients given by: 



/o 



— In In 

2 a + 5 



const. 



(28) 



/2 = _ (^±^ - I (a' - 5')(cos2cc>i + cos2c^2) -(29) 
4 o 

These expressions have been obtained in the limit d/Ac — > 
which, however, is attained very slowly. In this limit the 
free energy difference T{d^ uji) — T{d 00, Wi) is actually 
ill-defined and therefore the effective colloidal interaction 
is only meaningful for a finite capillary length Ac, simi- 
lar to a free two-dimensional interface the width of which 
is determined by the capillary wave fluctuations and di- 
verges logarithmically ~ In Ac- The leading terms of the 
fluctuation-induced free energy between spheres or disks 
(where the contact lines are circles of radius tq) have al- 
ready been calculated in Refs. [I91I20] and correspond to 
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the orientation-independent terms in Eqs. ([28|) and (j29p 
with 2ro = a + b. Note that for stretched eUipsoids with 
a > 2b the next-to-leading order term /2 in the free en- 
ergy expansion may become repulsive when the ellipsoids 
are aligned side by side (wi = UJ2 = 7''/2). 

In the opposite limit of small surface-to-surface dis- 
tance ho{d,toi,uj2) the fluctuation force can be calculated 
by using the well-known result for the fluctuation force 
per length f2d{h) = —k'QTTT/{12h'^) between two lines 
a distance h apart together with the Derjaguin (or 
proximity) approximation: 

^ 71-2 I 2 , . , , 

k^T 24 y /iQ (i?^ + i?2 ) 

Here, i?i and R2 are the radii of curvature at those points 
on the contact line of ellipsoid 1 and 2, respectively, whose 
distance is the minimal surface-to-surface distance /lo (see 
Fig.[5](b)). Thus, the fluctuation-induced interaction en- 
ergy between the ellipsoids diverges upon approach, simi- 
larly to the van-der-Waals attraction. 

For intermediate distances d the partition function must 
be evaluated numerically. In Eq. ([26|) the integral over the 
auxiliary fields can be carried out because they appear 
only quadratically in the exponent. The resulting deter- 
minant is divergent and requires regularisation. However, 
the derivative of its logarithm with respect to d (corre- 
sponding to minus the force Ffiuc in direction of the dis- 
tance vector between the centers of the ellipsoids) is finite 
and convergent in a numerical analysis (see Ref. [20j for 
further details). It turns out that the fluctuation force is 
attractive for all distances and orientations which were 
analyzed. This is already suggested by the close-distance 

regime (where i^fluc cx^ ~l/^o^^ always attractive) and 
the long-distance regime (where i^fiuc is dominated by the 
likewise attractive, in-plane isotropic term —dfo/dd = 
-l/[2dln{d/ro)], see Eq. ([28]) ). In order to exemplify the 
effect of in-plane anisotropy on the fluctuation force, the 
results for the force with the asymptotically leading, iso- 
tropic term subtracted (i^sub = -Ffluc + k^T dfo/dd) arc 
shown in Fig.[7]for ellipsoids with aspect ratio a/b ~ 6 and 
for the three configurations (a) tip-to-tip {uji = uj2 = 0°), 
(b) side-to-tip (wi = 90°, 102 = 0°) and (c) side-by-side 
(wi = UJ2 = 90°). For all configurations, for large d the 
approach to the aymptotic result given by —d{f2/d'^)/dd 
is fairly slow. For the configurations (a) and (b) the sub- 
tracted force Fsuh remains attractive for all distances and 
there is a smooth crossover from the longe-distance to 
the close-distance regime while for the side-by-side con- 
figuration (c) there is a sign change from the attractive 
close-distance regime (open circles) to the repulsive long- 
distance regime (full circles), in accordance with Eq. (^5]) . 

4 Summary and conclusions 

In this work we have analyzed the interface-mediated in- 
teractions which arise between ellipsoidal particles trapped 
at a fluid interface. 



Firstly, ellipsoids cause static interface deformations if 
they are partially wetting and their contact angle is differ- 
ent from 90°. These static deformations lead to orientation- 
dependent capillary interactions between the particles. The 
full solution to this capillary problem requires the solution 
of a nonlinear differential equation together with Young's 
condition on the boundary, the three-phase contact lines 
whose locations are a priori unknown. It is possible to 
analyze the interface deformation and the ensuing capil- 
lary potential in a perturbativc treatment, valid for small 
deformations of the interface, which leads to the standard 
problem of a linear differential equation with a local condi- 
tion on a given, fixed boundary. For small to intermediate 
distances between the ellipsoids we find considerable devi- 
ations from the well-known quadrupole interaction which 
is valid for asymptotically large distances. As a perspec- 
tive for future work, the developed algorithm allows a fast 
determination of the deformation and the potential also 
for large eccentricities of the particles and appears to be 
potentially useful for application in computer simulations 
of the aggregation process in ellipsoidal monolayers. 

Secondly, thermally excited capillary wave cause fluc- 
tuation-induced interactions between the ellipsoids. For 
the specific case of a pinned contact line we find that 
anisotropic effects in the fluctuation force arise only for 
subleading terms in an asymptotic expansion. It diverges 
for ellipsoids coming close to contact. However, due to its 
small scale the fluctuation force appears to be relevant 
experimentally only if the static capillary interactions are 
greatly reduced, e.g., if the ellipsoids are of nanometer size 
or the contact angle is close to 90°. 

Acknowledgment: E. Noruzifar and M. Oettel thank 
the German Science Foundation for financial support 
through the Collaborative Research Centre (SFB-TR6) 
"Colloids in External Fields" . 



A Contact line contributions to the free 
energy 

In this appendix we determine the coefficients Rzz^ Rza 
and Raa of the functional Taylor expansion of the bound- 
ary free energy in Eq. ([^0)1 around the reference config- 
uration {vi^^ = v{ro{(pi)) — Ahi = 0, ai.rof — 0}. They are 
given by the second variation of — Y^i=i with re- 
spect to shifts of the contact line height Vi^^. or to changes 
in the orientation ai of the long axis of ellipsoid i and 
can be calculated separately for the two colloids. As in 
Sec. 12.11 for the case of ai = (for a single ellipsoid), the 
needed functional derivative S^J-'h,i/6v'^^., the derivative 
d^^h,i/dai^ and the mixed derivative -^{5T'h,i/ 5vi^^.) 
are determined by the derivatives of the boundaries of the 
surface integrals (given by the position of the three phase 
contact line, cf. Eq. ^) with respect to Ui and Vi^^, re- 
spectively. 

According to Eq. ([7]), ~ ^{cosOAAi i -f AA^^a],i), 
upon contact line shift and tilt the boundary free energy 
contains a contribution due to the change of the air-water 
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interface (projected onto z = 0): 

/•27r rro{4>i) 

^^proj,» = - / del), drr , (31) 

and a contribution due to the change of the colloid area 
exposed to fluid I: 

AAi,^ = / dcj)' I dx' I dy' (32) 

X S \ (J)' — arctan ^ | a/ g(x' , y') 
V ^ J 

+ d(t)' dx' / dy' 



S ^<j)' — vr — arctan ^/g{x', y') 



In Eq. (|32|) . the surface integral is transformed into one 
over the cartesian components x', y' of a body-fixed coor- 
dinate system with axes fixed to the main ellipsoid axes 
for computational convenience in taking the derivatives. 
The reference contact line is parametrized by x^]~\(j)'), 
y^}~\4>') (for (/)' < TT [(/)' > tt]), and the shifted contact 
line is parametrized by x'+I"' (0'), y'+'^l (0'). 

The calculation of the boundary free energy variation 
^b,i can be performed in two steps: (1) First, we tilt el- 
lipsoid i by an angle with a pinned contact line. Then 
AAi^i = holds, but AAp^a],i 7^ 0. (2) In a second step, 
the contact line is released to its final position Vi^^- . In this 



second step, both AAi^i and AA 



proj.z, 



contribute to Th 



The two steps have to be distinguished, since the orien- 
tational tilts change the surface measure of the ellipsoid. 
In order to avoid the calculation with the a^-dependent 
metric g = 1 -t- |Vzi^Eii(aj)P of tf^c ellipsoid surface, we 
calculate the second variation S^AAi^i in body-fixed co- 
ordinates of the colloid. The contribution 6'^AAp^oj,i from 
the change of the projected meniscus area is determined 
separately, here the steps (1) and (2) can be considered 
together. In body-fixed coordinates the metric is given by 
g'^l + {l-eYix'/z'f + {y'/z'r. 

With Zi^^ = Vi^^^ + u^ci{4'i) - /ii,rcf being the height 
of the contact line relative to the center of colloid i, the 
position of the three contact line in body-fixed coordinates 
after a shift w^^^. of the contact line and a rotation of the 
ellipsoid by an angle ai (with respect to the plane z = 0) 
is determined by the equations 



x' ^ z sin a + a;((/), a) cos a , 
z' — z cos a — x{4>^ a) sin a , 

y' 

parametrising the rotation, and the constraint 



r 1,2 12 /2 

t ~ b — z — y 



{l-e')x'' 







(33) 
(34) 
(35) 



(36) 



which ensures that the contact line is on the ellipsoid sur- 
face. 



Employing the relations given above we can calculate 
the contact line position on the rotated ellipsoid and, in 
particular, its partial derivatives. After some algebra, we 
finally arrive at the expressions 



Rzz{ 
Rza{ 

and 



Raa{4>i) — 
1 



52 sin^ e 



^0, 



^) = - 



cf(0.)'(l 
1 



cos^ (piY ' 

^0,rcf(4)^/''0,rof(</'i)^"' 



(37) 



1 — cos^ 



1 



1 



cos2 [(1 - e2) -f 6Vro,ref (0»)' 



(1 



(38) 



(39) 



1 — cos^ (f. 

cos^ 9 



ro,rol{(l>if COS^ ( 



4 set (1 + cos^ , 
1 — cos^ (hi 



2 (1 - e2 cos2 0,)3 
-(1 - e2)(62 _ zo^rcf (0»)') (1 + cos^ 0,(1 - 2e^ sin^ 0,)) 

+ 2&2 (1 - cos^ (t>^) (1 + cos^ 0J (1 + sir? 0,)) , 

for the coefficients of the functional Taylor expansion of 
JS. Here, zo,rcf((/>)^ = 9 -ro^,^f{4>f[ii\T? 0-f (l-e^) cos^ 0]. 

Note, that by minimizing Th with respect to ai we 
ensure torque balance of the colloid whereas mimimizing 
with respect to the height Ahi guarantees that the to- 
tal vertical force exerted by the meniscus on the contact 
line vanishes. The numerical calculation of the equilibrium 
meniscus between two ellipsoidal particles and the result- 
ing capillary interactions shows, however, that the effect of 
ai is rather small, leading to changes of the results ^1% 
as compared to the computations using the boundary con- 
tribution ([13]), in which ai is neglected. 
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Fig. 7. Results for the fluctuation force with the leading 
asymptotic term subtracted, -Fsub ~ Jfluc + ksT dfo/dd, for 
ellipsoids with aspect ratio a/b = 6 and for the three conflgura- 
tions (a) tip-to-tip {uii = lj2 — 0°), (b) side-to-tip (oji = 90°, 
UJ2 =0°) and (c) side-by-side (oji — ui2 ~ 90°). Numerical re- 
sults are shown by circles, the next-to-leading asymptotic term 
involving the coefflcient /2 (Eq. (|29|) ) is represented by a full 
line, and the Derjaguin approximation (derivative of Eq. (|30p 
with respect to ho) is given by a dashed line, respectively. The 
capillary length was chosen as Ac — 10^ b. 



